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ABSTRACT 

The origin of the recently reported systematic bias in the spectroscopic temperature of galaxy 
clusters is investigated using cosmological hydrodynamical simulations. We find that the local inho- 
mogeneities of the gas temperature and density after corrected for the global radial profiles, have 
nearly a universal distribution that resembles the log-normal function. Based on this log-normal ap- 
■ proximation for the fluctuations in the intra-cluster medium, we develop an analytical model that 

explains the bias in the spectroscopic temperature discovered recently. We conclude that the multi- 
phase nature of the intra-cluster medium not only from the radial profiles but also from the local 
inhomogeneities plays an essential role in producing the systematic bias. 

Subject headings: galaxies: clusters: general - X-rays: galaxies masses - cosmology: observations 

a' 

1. INTRODUCTION 

| Recent progress both in numerical simulations and observations has improved physical modeling of galaxy clus- 
. ters beyond a simple isothermal and spherical approximation for a variety of astrophysical and cosmological ap- 
00 ' plications; departure from isothermal distribution was discussed in the context of the S unyaev-Zel'dovich effect 
(jlnagaki. Suginohara fc Sutol 119951 : lYoshikawa. Itoh fc Sutol 119981 : lYoshikawa fc Sutol I1999D. an empirical 0- -model 



X 



profile has been replaced by that based on the NFW dark matter density profile ((Navarro. Frenk. fc White! 119971 ; 
iMakino. Sasaki fc Sutol 1998; Suto. Sasaki. & Makino 19 9jj), non-spher i cal effects of dark halos have been consid- 
ered fairly extensively (|Lee fc Shandarinl Il998t ISheth fc Tormenl 1 19991 Ping fc Sutol l2002t iLee fc Sutol 12001 12004 
■ iKasun fc Evrardl l2005h . and the physical model for the origin of the triaxial density profile has been proposed 
Q ■ (|Lee. Jing fc Sutdl2005h . 



Despite an extensive list of the previous studies, no physical model has been proposed for the statistical nature of 
~ q ■ underlying inhomogeneities in the intra-cluster medium (ICM, hereafter). Given the high spatial resolutions achieved 
both in observations and simulations, such a modeling should play a vital role in improving our understanding of 
O | galaxy clusters, which we will attempt to do in this paper. 

+-> ■ Temperature of the ICM is one of the most important quantities that characterize the cluster. In X-ray observations, 
the spectroscopic temperature, T spcc , is estimated by fitting the thermal continuum and the emission lines of the 
. . , spectrum. In the presence of inhomogeneities in the ICM, the temperature so measured is inevitably an averaged 
quantity over a finite sky area and the line-of-sight. It has been conventionally assumed that T spec is approximately 
equal to the emission-weighted temperature: 

_ Jn 2 A(T)TdV 

° w - Jn 2 A(T)dV ' ( ' 

where n is the gas number density, T is the gas temperature, and A is the cooling function. Ma zzotta et al.l {2004), 
however, have pointed out that T spec is systematically lower than T cw . The authors have proposed an alternative 
definition for the average, spectroscopic-like temperature, as 

= Jn 2 T a - 1 / 2 dV 

Tsl = Jn 2 T°-VW (2) 

They find that T B \ with a = 0.75 reproduces T spec within a few per cent for simulated clusters hotter than a few keV, 
assuming Chandra or XMM-Newton detector response functions. Ra sia etldl (|2005l ) performed a more systematic 
study of the relation between T nw and T s \ using a sample of clusters from SPH simulations and concluded that 
T B \ ~ 0.7T CW . IVikhlininl ((2006) provided a useful numerical routine to compute T s \ down to ICM temperatures of 
~ 0.5 keV with an arbitrary metallicity. It should be noted that T ew is not directly observable, although it is easily 
obtained from simulations. 
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The above bias in the cluster temperature should be properly ta ken into account w hen confronting observational 
data with theory, for example, in cosmological st udies. As noted bv iRasia et al.l (|2005h . it can result in the offset in 
the mass-temperature relation of galaxy clusters. iShimizu et al.l (|2006l ) have studied its impact on the estimation of 
the mass fluctuation amplitude at 8/i _1 Mpc, erg. The authors perform the statistical analysis using the latest X-ray 
cluster sample and find that erg ~ 0.76 ± 0.01 + 0.50(f — «m), where ajf = T spec /T ew . The systematic difference of 
Tspec ~ 0.7T ew can thus shift erg by ~ 0.15. 

In this paper, we aim to explore the origin of the bias in the spectroscopic temperature by studying in detail the 
nature of inhomogeneities in the ICM. We investigate both the large-scale gradient and the small-scale variations 
of the gas density and temperature based on the cosmological hydrodynamical simulations. Having found that the 
small-scale density and temperature fluctuations approximately follow the log-normal distributions, we construct an 
analytical model for the local ICM inhomogeneities that can simultaneously explain the systematic bias. 

The plan of this paper is as follows. In §2, we describe our simulation data, construct the mock spectra and 
compare quantitatively the sp ectroscopic temperatur e with the emission- weighted temperature and the spectroscopic- 
like temperature suggested bv lMazzotta et al] (|2004T ). In §3, we propose an analytical model for the inhomogeneities 
in the ICM. In §4, we test our model against the result of the simulation. Finally, we summarize our conclusions in 
§5. Throughout the paper, temperatures are measured in units of keV. 

2. THE BIAS IN THE SPECTROSCOPIC TEMPERATURE 
2.1. Cosmological Hydrodynamical Simulation 

The results presented in this paper have been obtained by using the fina l outp ut of the Smoothing Particle Hydro- 
dynamic (SPH) sim ulation of the l ocal u niverse performed by iDolag et al.l ((2005). The initial conditions were similar 
to those adopted by iMathis et~a l. (2002) in their study (based on a pure N-body simulation) of structure formation 
in the local universe. The simulation assumes the spatially-flat A cold dark matter (A CDM) universe with a present 
matter density parameter Qom = 0.3, a dimensionless Hubble parameter h — 7?o/100 km s _1 Mpc -1 = 0.7, an rms 
density fluctuation amplitude erg = 0.9 and a baryon density parameter fib = 0.04. Both the number of dark matter 
and SPH particles are ~ 50 million within the high-resolution sphere of radius ~ 110 Mpc which is embedded in a 
periodic box of ~ 343 Mpc on a side, filled with nearly 7 million low-resolution dark matter particles. The simulation 
is designed to reproduce the matter distribution of the local Universe by adopti ng the initial conditions based on the 
IRAS galaxy distribution smoothed on a scale of 7 Mpc (seelMathis et al.l l2002. for detail). 

The run has been carried out with GADGET-2 llSpringel 12005^ 7 a new version of the parallel Tree-SPH simulation 
code GADGET (|Springel et al.| [2001). The code uses an entropy-conserving formulation of SPH (|Springel &; Hernquistj 
2002), and allows a treatment of radiative cooling, heating by a UV background, and star formation and feedback 
processes. The latter is ba sed on a sub- resolution model for the multiphase structure of the interstellar medium 
(|Springel fc Hernquistj |2003): in short, each SPH particle is assumed to represent a two-phase fluid consisting of cold 
clouds and ambient hot gas. 

Th e code also follows the pattern of metal production from the past history of cosmic star formation (jTornatore et al.l 
2004). This is done by computing the contributions from both Type-II and Type-la supernovae and energy feedback 
and metals are released gradually in time, accordingly to the appropriate lifetimes of the different stellar populations. 
This treatment also includes in a self-con sistent way the dependence of the gas cooling on the local metallicity. The 
feedback scheme assumes a Salpeter IMF (|Salpeterlll955| ) and its parameters have been fixed to get a wind velocity of 
ps 480 km s _1 . In a typical massive cluster the SNe (II and la) add to the ICM as feedback w 2 keV per particle in an 
Hubble time (assuming a cosmological mixture of H and He); w 25 per cent of this energy goes into winds. A more 
detailed discussion of cluster properties and metal d istribution within the IC M as resulting in simulations including the 
metal enrichment feedback scheme can be found in iTornatore et all (|2004f ). The simulation provides the metallicities 
of the six different species for each SPH particle. Given the fact that the major question that we addressed is not the 
accurate estimate of T spcc or T s \, but the systematic difference between the two, we decided to avoid the unnecessary 
complication and simply to assume the constant metallicity. Therefore we adopt a constant metallicity of 0.3 Zq in 
constructing mock spectra below, and the MEKAL (not VMEKAL) model for the spectral fitting. 

The gravitational force resolution (i.e. the comoving softening length) of the simulations has been fixed to be 14 kpc 
(Plummer-equivalent) , which is comparable to the inter-particle separation reached by the SPH particles in the dense 
centers of our simulated galaxy clusters. 

2.2. Mock Spectra of Simulated Clusters 

Among the most massive clusters formed within the simulation we extracted six mock galaxy clusters, contrived to 
resemble A3627, Hydra, Perseus, Virgo, Coma, and Centaurus, respectively. Table Q] lists the observed and simulated 
values of the total mass and the radius of these clusters. In order to specify the degree of the bias in our simulated 
clusters, we create the mock spectra and compute T spcc in the following manner. 

First, we extract a 3/i _1 Mpc cubic region around the center of a simulated cluster and divide it into 256 3 cells so 
that the size of each cell is approximately equal to the gravitational softening length mentioned above. The center 
of each cluster is assigned so that the center of a sphere with radius l/i _1 Mpc equals to the center of mass of dark 
matter and baryon within the sphere. 
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The gas density and temperature of each mesh point (labeled by I) are calculated using the SPH particles as 

Ngas 

Pi= mW(\ri-ri\,hi), (3) 
i=i 

Ngas _, 

Tz= V -^iwQn-nlhi), (4) 
,-=i Pl 

where rj is the position of the mesh point, IF denotes the smoothing kernel, and mj, r^, TJ, and are the mass 
associated with the hot phase, position, smoothing length, temperature, and density associated with the hot gas phase 
of the i-th SPH particle, respectively. We adopt the smoothing kernel: 

1 ( 1 - (3/2)u 2 + (3/4)u 3 if < u < 1 
W{\ri-r i \,h i ) = — s \{2-uf /A if 1 < u < 2 (5) 

ntl i [0 otherwise, 

where u = |rj — ri\/hi. 

It should be noted that the current implementation of the SPH simulation results in a small fraction of SPH particles 
that have unphysical temperatures and densities. This is shown in the temperature - density scatter plot of Figure 
[TJ The red points correspond to SPH particles that should be sufficiently cooled, but not here because of the limited 
resolution of the simulation. Thus if they satisfy the Jean criterion, t hey should be re g arded as simply cold clumps 
without retaining the hot gas nature; see Figure 1 and section 2.1 of lYoshikaw a et al.l ([2001). In contrast, the blue 
points represent the SPH particles that have experienced the cooling catastrophe, and have significantly high cold gas 
fraction (larger than 10 percent). In either case, they are not supposed to contribute the X-ray emission. Thus we 
remove their spuriou s contribution to the X-ray emission and the temperature estimate of ICM. Specifically we follow 
iBorgani et al.l (|2004l ). and exclude particles (red points) with 7$ < 3 x 10 4 K and pi > 500p c 51b , where p c is the critical 
density, and particles (blue points) with more than ten percent mass fraction of the cold phase. While the total mass 
of the excluded particles is very small (~ 1%), they occupy a specific region in the p-T plane and leave some spurious 
signal due to high density, in particular for blue points. 

Second, we compute the photon flux f(E) from the mesh points within the radius r2oo from the cluster center as 

f(E)dE ocexp (-a gal {E)N H ) V - P } , 4 ( P om {T u Z, E(l + z cl )) dE(l + z cl ), (6) 

4tt 1 + z c i) 4 \m*J 

where z c \ denotes the redshift of the simulated cluster, Z is the metallicity (we adopt 0.3 Zq), X is the hydrogen 
mass fraction, m p is the proton mass and P em (Tj,Z,E) is the emissivity assuming collisional ionization equilibrium. 
We calculate P em (Ti,Z,E) using SPEX 2.0. The term exp (— o-gaiNn) re presents the galactic extinction- Nh is the 
column density of hydrogen and a ga i{E) is the absorption cross section of iMorrison fc McCammonl |T983). Since we 
are interested in the effect due to the spectrum distortion, not statistical error, we adopt a long exposure time as 
the total photon counts TV = J^Zq° 5 Ef(E) dE ~ 500,000. In this paper, we consider mock observations using 
Chandra and XMM-Newton, thus we neglect a peculiar velocity of the cluster and a turbulent velocity in ICM because 
of insufficient energy resolution of Chandra ACIS-S3 and XMM-Newton MOS1 detector. 

Finally, the mock observed spectra are created by XSPEC version 12.0. We consider three cases for the detector 
response corresponding to 1) perfect response, 2) Chandra ACIS-S3, and 3) XMM-Newton MOS1. In the first case, 
we also assume no galactic extinction (Nh = 0) and refer to it as an "IDEAL" case. In the second and third cases, we 
adopt an observed value to Nh listed in Table Q] and redistribute the photon counts of the detector channel according 
to RMF (redistribution matrix file) of ACIS-S3 and MOS1 using the rejection method. 

Figure [2] illustrates the mock spectra of "Virgo" and "Perseus" using RMF of ACIS-S3. Unless stated otherwise, 
we fit the spectra by an absorbed single-temperature MEKAL model in the energy band 0.5-10.0 keV. We define the 
spectroscopic temperature, T spoc , as the best- fit temperature provided by this procedure. Since the spatial resolution 
of the current simulations is not sufficient to fully resolve the cooling central regions, a single-temperature model yields 
a reasonable fit to the mock spectra. For comparison, we also plot the spectra for a single temperature corresponding 
to the "emission weighted" value of the mesh points within r2oo : 

We calculate the cooling function A(T) using SPEX 2.0 assuming collisional ionization equilibrium, the energy range 
of 0.5-10.0 keV, and the metallicity 0.3Z Q . The difference between T spoc and T c s ™ ,m is clearly distinguishable on the 
spectral basis in the current detectors. 

Figure [3] shows the relation between T spoc and T ew for our sample of simulated clusters. It is well represented by a 
linear relation T spoc = kT cw + I with the fitted values of k = 0.84, Z = 0.34 (IDEAL), k = 0.84, 1 = 0.36 (ACIS) and 
k = 0.85, 1 = 0.31 (MOS), respectively. In the range of temperatures corresponding to rich clusters, the spectroscopic 
temperature T spec is systematically lower than T cw by 10 — 20 %. 
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We note that the above bias should depend on the energy band in which T spoc is evaluated. In order to demonstrate 
it quantitatively, we also list in Table Q] the fitted values of T spcc from the 0.1-2.4 keV and 2.0-10.0 keV data, respec- 
tively. Because the exponential tail of the thermal bremsstrahlung spectrum from hotter components has negligible 
contribution in the softer band, the bias tends to increase and decrease in the softer and harder bands, respectively. 

2.3. Spectroscopic- Like Temperature 

In order to better approximate T spe c, Mazzott a et al.l (|2004h proposed a " spectroscopic-like temperature"; they 
found that equation (|2|) with a = 0.75 reproduces T spe c in the 0.5-10.0 keV band within a few percent. Throughout 
this paper, we adopt a = 0.75 when we estimate the spectroscopic-like temperature quantitatively. In Figure [31 we 
also plot this quantity computed from the mesh points within r2oo- 

2 T 0.25 

rrisim,m L-iI^r wa "I I lo\ 

si — ^ „2^-0.75 • \°> 

^ier 200 Pi 1 1 

As indicated in the bottom panel, T^ m ' m reproduces T spec within 6 % for all the simulated clusters in our sample. 
Given this agreement, we hereafter use T^ m ' m to represent T spec , and express the bias in the spectroscopic temperature 

by 

_ rrisim,m ;rrisim,m /q\ 
fCsim — si / J cw ■ V y 7 

Table [T] provides K S i m (mesh-wise) for the six simulated clusters. The range of K s i m (mesh-wise) is ap proximately 
0.8-0.9. While 

K sim is systematically lower than unity, the value is somewhat higher than the results of iRasia et al.l 
(|2005f ). T s i ~ 0.7T CW . This is likely due to the different physics incorporated in the simulations and the difference in 
how T s i and T cw are compu ted from the simu lation outputs. The major difference of the physics is the amplitude of 
the wind velocity employed: IRasia et alj (|2005l ) used have a feedback with weaker wind of 340km s _1 , while our current 
simulations adopt a higher value of 4 80km s -1 . Because weaker wind cannot remove small cold blobs effectively, the 
value of T s i/T ew of IRasia et al.l (|2005l ) is expected to be larger. 

To show the difference of the temperature computation scheme explicitly, we also list in Table Q] the values of 
K s i m = T s S j lm ' p /T s ™' p (particle- wise) computed in the "particle- wise" definitions used in Ras ia et al.l ((2005): 

T sim,p _ Sigr 200 mjPik(Ti)Ti 



and 



rr^Sim,p 



£ ieraoo ™^ - 2 '' 



(see also lBorgani et al1l2004l) . In practice, the emission- weighted and spectroscopic-like temperatures defined in equa- 
tion (fT0|) and equation (fTTj) are sensitive to a small number of cold (and dense) SPH p articles present, w hile their 
contribution is negligible in the mesh- wise definitions, equation ((7]) and equation ((§]). As in Ra sia et"aT1 (|2005f L we have 
removed the SPH particles below a threshold temperature Tu m = 0.5 keV to compute K s ; m (particle- wise) . Table [1] 
indicates that K S j m (particle- wise) tends to be systematically smaller than K S j m (mesh- wise). Even adopting Tn m = 0.01 
keV makes K s i m (particle- wise) smaller only by a few percent. 

Given the limit of the particle-wise definitions mentioned above, we use the mesh-wise definitions of the emission- 
weighted temperature (T c s ™' m ) and the spectroscopic-like temperature (T^ m ' m ) given in equation (J7J) and (JSJ), respec- 
tively, in the following sections. 

3. ORIGIN OF THE BIAS IN THE SPECTROSCOPIC TEMPERATURE 
3.1. Radial Profile and Log-normal Distribution of Temperature and Density 

Having quantified the bias in the temperature of simulated clusters, we investigate its physical origin in greater 
detail. Since clusters in general exhibit inhomogeneities over various scales, we begin with segregating the large-scale 
gradient and the small-scale fluctuations of the gas density and temperature. 

For the large-scale gradient, we use the radially averaged profile of the gas temperature and density shown in FigurelH 
We divide the simulated clusters into spherical shells with a width of 67ft. -1 kpc and calculate the average temperature 
T(r) and density n(r) in each shell. The density profile n(r) is fitted to the conventional beta model given by 

n(r) = n 



1 + (r/r c ) 



(12) 



where no is the central density, r c is the core radius, and f3 is the beta index. We adopt for the temperature profile 
T(r) the polytropic form of 

T(r) = T [n(r)/n o r~\ (13) 

where To is the temperature at r = 0, and 7 is the poly trope index. The simulated profiles show reasonable agreement 
with the above models. The best-fit values of j3 and 7 are listed in Table [1] The range of 7 is approximately 1.1-1.2. 
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In addition to their radial gradients, the gas density and temperature have small-scale fluctuations. Figure [5] 
illustrates the distributions of the gas density and temperature in each radial shell normalized by their averaged 
quantities, (T) and (n), respectively. Despite some variations among different shells, we find a striking similarity in the 
overall shape of the distributions. They approximately follow the log-normal distribution given by 

,2- 



X 



1 



27TCr L N,x 



exp 



(log 5 X 



'LN 



■/2)' 



dS x 



(14) 



where 8 X = x/{x) and x denotes T or n (St = T/(T),8 n = n/(n)). For simplicity, we neglect the variations among 
different shells and fit the distribution for the whole cluster within r 2 oo (solid line) by the above equation (dashed 
line). The best-fit values of ctln.t and ctln,™ are listed in Table [T] 

The small-scale fluctuations mentioned above are not likely an artifact of the SPH scheme. We have applied the 
similar analysis to the data of grid-based simulations (D.Ryu, private communication) and obtained essentially the 
same results. Thus the log-normal nature of the fluctuations is physical, rather than numerical. 

3.2. Analytical Model 

Based on the distributions of gas density and temperature described in §3.1, we develop an analytical model to 
describe the contributions of the radial profile (RP) and the local inhomogeneities (LI) to the bias in the spectroscopic 
temperature. 

To describe the emission-weighted and spectroscopic-like temperatures in simple forms, let us define a quantity: 



A a = J n(r) 2 T(r) a dr, 



r 2 dr J dnn 2 {r)T a (r), (15) 

where the second line is for the spherical coordinate and R denotes the maximum radius considered here (we adopt 
R = ?'2oo)- Using this quantity, we can write down T s \ via equation @ as 



1 s\ 



del 



A a -i /2 /A a _ 3/2 . (16) 

When the temperature is higher than ~ 3 keV, the cooling function is dominated by the thermal bremsstrahlung; 
A cx s/T. We find that substitution \ff for A(T) yields only ~ 1% difference for the simulated clusters. In the present 
model, we adopt for simplicity the thermal bremsstrahlung cooling function, A oc VT, and then equation ([1]) reduces 
to 

T™ dcl = A 3/2 /A 1/2 . (17) 
When evaluating equation (|15[) . we replace the spatial average with the ensemble average: 



dnn 2 (r)T a (r)=4n{(n)(r)Y[(T}(r)] a / 5 z J^P{5 n ,5 T -r)d8 n d5 T , 



(18) 



where P(S n , St', t) is a joint probability density function at r. Assuming further that the temperature inhomogeneity 
is uncorrelated with that of density, i.e. P(6 n , Sr',r) = P(5 n ; r) P(St; r), we obtain 



A, 



An 



r%n)(r)Y[(T){rWdr d 2 n P(S n ;r)dS. 



8%P{8 T ;r)d5 1 



(19) 



n 



For the log-normal distribution of the temperature and the density fluctuations, the average quantities are expressed 
as 

~a(a - l)cr LNiT (r) 2 



8^P(5t] r)d8r = exp 



SlP(S n ; r)d8 n = exp [cln, n{ r ) z 



(20) 
(21) 



If <7ln, t and <TLN,n are independent of the radius (ctln, t(p) — 0ln,t , cln, n( r ) — cln, n), equation (JTUJ) reduces to 



A a = exp(of jN ) exp 



a(a — 1) 



2 

CT LN,T 



X 47T 



r 2 [(n)(r)] 2 [<T)(r)]^r. 



Using the above results, 7 1 ™ " 301 and T c ™ odel are expressed as 

T s r dcl = A I _ 1/2 A4 a _ 3 /2 = ^exp 



a 2 / ^LN.^ 



rpll 



■A 3/2 /A 1/2 = T^exp [ -cr£ NjT 



(22) 

(23) 
(24) 



6 



Kawahara et al. 



where T RP and T e ^f are denned as 



,. tfr 2 [(n)(r)] 2 [{T)(r)](°-Wdr 



t hp = ^ ivv,, u /v ^ » (25) 

/„ r 2 [(n)(r)] 2 [(T)(r)]( a - 3 / 2 )dr 

t rp == / V [(»)(r)] 2 [(r)(r)]3/ 2 dr 
/ V [<n)(r)] 2 [(T)(r)]Va dr ' 

As expected, r™ odcl and X'™ do1 reduce to T RP and Tj^ 3 in the absence of local inhomogeneities (<tln,t = 0). Note 
that equations (f23|) and (|2~t|) are independent of ctln, n- This holds true as long as the density distribution P(S n ) is 
independent of r. 
The ratio of T s f odel and T™ dcl is now written as 

„ — fpmodel /remodel RP LI I07\ 

"-model — 1 si / cw 5 

where k rp and k li denote the bias due to the radial profile and the local inhomogeneities, 

k rp _ r RP /T RP, (28) 



and 



K LI = exp [(a - 2) cr£ N)T ] , (29) 

respectively. Figure [6] shows k li for a fiducial value of a = 0.75 as a function of oln. t- The range of <tln, t of simulated 
clusters, 0.1 < <tln,t < 0.3, (Table [J) corresponds to 0.99 > n LI > 0.89. 

In the case of the beta model density profile (eq.[l2]) and the polytropic temperature profile (eq.QjB]), T RP and T™ 1 
are expressed as 

K ™=T*?/T™, (30) 
rpRp _ 2^1 (3/2, 3g[l + ( 7 - l)(a - l/2)/2]; 5/2; ~P 2 /r 2 c ) 

sl 2 ^(3/2, 3/3[l+ ( 7 -l)(a-3/2)/2]; 5/2; -i? 2 /r c 2 ) °' 

RP _ 2 Fi(3/2, 3/?[l + 3(7 - l)/4]; 5/2; -R 2 /r 2 ) 

ow 2^(3/2,3/311 + ( 7 - l)/4]; 5/2; -i? 2 /r 2 ) °' 1 ' 

where 2-Fi(a, /3; 7; C) is th e hyper geometric function. 

Figure [7] shows k rp as a function of f3 for various choices of 7 and r c /r2oo- Given th at a number of observed clusters 
exhibit a cool core, we also plot the case with the temperature profile of the form ([Allen. Schmidt fc Fabiar] 120011 : 
iKaastra et aLll2004l) : 

T(r) =T l + (T h T,) -^f-, (33) 

with (Th — Ti)/T = 1.5 and fi = 2. For the range of parameters considered here, k rp exceeds 0.9. This implies that 
the bias in the spectroscopic temperature is not fully accounted for by the global temperature and density gradients 
alone; local inhomogeneities should also make an important contribution to the bias. 

4. COMPARISON WITH SIMULATED CLUSTERS 

We now examine the extent to which the analytical model described in the previous section explains the bias in the 
spectroscopic temperature. The departure in the radial density and temperature distributions from the beta model 
and the polytropic model results in up to 7 % errors in the values of and T S RP . Since our model can be applied to 
arbitrary (n) (r) and (T) (r) , we hereafter use for these quantities the radially averaged values calculated directly from 
the simulation data. We combine them with cr LN , T in Table ED to obtain T™ odo1 (eq.[23]) and T c ™ odcl (eq.^J). 

Figure E compares T™ odcl and T" odcl against T^ m ' m and T^ m (eq.0 and eq.0), respectively. For all clusters 
except "Perseus" , the model reproduces within 10 percent accuracy the temperatures averaged over all the mesh points 
of the simulated clusters. 

Given the above agreement, we further plot Kmodei against K s i m in Figure [9] The difference between K S i m and Kmodei 
is kept within ~ 10 % in all the cases. Considering the simplicity of our current model, the agreement is remarkable. 
In all the clusters, both k rp and k li are greater than K sml , indicating that their combination is in fact responsible for 
the major part of the bias in the spectroscopic temperature. 

In §3.2, we assumed that n and T are uncorrelated, i.e., P(£„,<$t) = P{S n )P(ST)- We examine this assumption 
in more detail. We pick up two clusters, "Hydra" and "Perseus", which show the best and the worst agreement, 
respectively, between K mo dei and K sml . Figure flOl shows the contours of the joint distribution of St = T(r) / '(T)(r) and 
S n = n(r)/(n)(r) for all the mesh points within r2oo for these clusters, together with that expected from the model 
assuming P(<j„,&r) = Pln(Sti)Pln(St)- For the log-normal distributions, Pln(^ti) and Pln{St), we have used the 
fits shown as the dashed line in Figure [5] The joint distribution agrees well with the model for both cases, while the 
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deviation is somewhat larger in "Perseus" for which the model gives poorer fits to the underlying temperature and 
density distribution in Figure [5] If the cluster is spherically symmetric and the ICM is in hydrostatic equilibrium, we 
expect that n are correlated with T as S u St = 1- We do not find such correlations in Figure [TU1 

Although the effects of the correlation is hard to model in general, we can show that it does not change the value of 
Kmodei = T s I f odol /T c ™ odcl as long as the joint probability density function follows the bivariate log-normal distribution: 



(1-p' 2 )- 1 / 2 
^bln (yn j S T ) dS n dS T = - exp 

^7raLN,nO"LN,T 



A 2 - Ip'AB + B 2 



2(1 -p 



dS n <1St 



(34) 



where p' = log [p(cxp cr^ N „ - l) 1 / 2 (exp^ NiT - 1) 1/2 + l]/cr LN , tctln, n , A = log(<5„) + cr£ N) „/2, B = log(<5 T ) + 
(t£ n t/2, and p is the correlation coefficient between n and T. Adopting p — yields pBLN{S n , St) = Fln^h) Pln(St)- 
The marginal probability density function of density J gJ^Pbln^ti! St) and that of temperature J d5 n PBLN(S n , St) 

^modcl Trao' 
si ' cw 

r 



are equal to P L n(^) and Pm(S T ), respectively. Using PblnO^, <5 T ), we obtain T™ odel , T™ odel as 



remodel _ rpRP 
1 si — 1 si ex P 



2 



^ °LN,T + 2/9 / CTLN,TO'LN, 



(35) 



1 



T c r dcl = exp i^-a 2 Nt T + 2pV LN , t ^ln, „J . (36) 

Although both T™ dcl and T™ do1 increase with the correlation coefficient, K modol = T^ odel /T™ del remains the same 
as that given by equation (f2"9)) . 

5. SUMMARY AND CONCLUSIONS 

We have expl o red th e origin of the bias in the spectroscopic temperature of simulated galaxy clusters discovered by 
Mazz otta et al.l ([2004D . Using the independent simulations data, we have constructed mock spectra of clusters, and 
confirmed their results; the spectroscopic temperature is systematically lower than the emission-weighted temperature 
by 10-20% and that the spectroscopic-like temperature defined by equation ([2]) approximates the spectroscopic tem- 
perature to better than ~ 6%. In doing so, we have found that the multi-phase nature of the intra-cluster medium is 
ascribed to the two major contributions, the radial density and temperature gradients and the local inhomogeneities 
around the profiles. More importantly, we have shown for the first time that the probability distribution functions of 
the local inhomogeneities approximately follow the log-normal distribution. Based on a simple analytical model, we 
have exhibited that not only the radial profiles but also the local inhomogeneities are largely responsible for the above 
mentioned bias of cluster temperatures. 

We would like to note that the log-normal proba bility dis t ribution functions for den s ity fields show up in a 
variety of astrophysical/co s mological problems (e.g., iHubble I Il934t iColes &: Jones! 119911 : iWada fc Normanl 120011 : 
iKavo. Taruva fc Sutol l200ll iTaruva et all l2002f ). While it is not clear if they share any simple physical principle 
behind, it is interesting to attempt to look for the possible underlying connection. 

In this paper, we have focused on the difference between spectroscopic (or spectroscopic-like) and emission-weighted 
temperatures, which has the closest relevance to the X-ray spectral analysis. Another useful quantity is the mass- 
weighted temperature defined by 

/ nTdV 

mw = ~J^dV~' (37) 

This is related to the cluster mass more directly (e.g..[Mathicsc n~fc Evrarc3l2001l : iNagai. Vikhlinin. fc Kravt sov 2006). 
The mass- weighted temperature is highly sensitive to the radial density and temperature profiles, while it is little 
affected by the local inhom ogeneities. Though challenging, it will be "observable" either by high-resolution X-ray 
spectroscopic observations (Na gai. Vikhlinin. fc Kravtsovl l2006f) or by a combination of the lower resolution X -ray 
spectroscopy and the Sunyaev-Zel'dovich imaging observations (|Komatsu et al.lll999L l2001: Kit avama et~aTl l2004) . We 
will discuss usefulness of this quantity and the implications for the future observations in the next paper. 
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TABLE 1 

Properties of the six simulated clusters and observed clusters. 



Simulation 




A3627 


Hydra 


Perseus 


Virgo 


Coma 


Centaurus 


M 2 oo[W 14 h- 1 M e ] 


2.2 


1.8 


6.7 


3.1 


4.3 


2.5 


r 2 oo [h' 1 Mpc] 


1.1 


1.0 


1.6 


1.2 


1.4 


1.1 


P 


0.69 


0.72 


0.65 


0.59 


0.73 


0.73 


7 


1.21 


1.26 


1.09 


1.12 


1.19 


1.14 


T ew [keV] 


4.2 


4.1 


5.7 


3.7 


6.7 


4.2 


T sl [keV] 


4.0 


.3.8 


4.7 


3.2 


5.9 


3.9 


T Bpcc (IDEAL) [keV] (0.5-10.0 keV) 


4.1 


.3.9 


5.0 


3.3 


6.1 


3.9 


T spcc (ACIS) [keV] (0.5-10.0 keV) 


4.1 


.3.9 


5.0 


3.3 


6.0 


4.0 


Tspcc (MOS) [keV] (0.5-10.0 keV) 


4.0 


3.8 


5.0 


3.3 


6.0 


4.0 


T spcc (IDEAL) [keV] (0.1-2.4 keV) 


1.0 


.3.7 


4.7 


3.1 


5.9 


3.9 


T S pcc (IDEAL) [keV] (2.0-10.0 keV) 


4.1 


4.0 


5.4 


3.6 


6.5 


4.0 


ftsim (mesh-wise) 


0.95 


0.92 


0.84 


0.86 


0.88 


0.95 


K sim (particle-wise) 


0.88 


0.89 


0.70 


0.75 


0.84 


0.86 




0.97 


0.94 


0.98 


0.94 


0.96 


0.97 


CT LN, T 


0.159 


0.180 


0.316 


0.286 


0.159 


0.178 


CLN, n 


0.240 


0.180 


0.518 


0.446 


0.434 


0.239 



Observation 



A3627 



4 -°-0.58 
*1.26 
5 62+ - 12 
21.7 



Hydra 



Perseus 



1 qn+ l, ' ;i8 
1-MU_ 33 

1.22 

W.20 



9-08+f.52 
2.05 

f0.06 



Virgo 



Coma 



Centaurus Ref 



TUT23 - 



M 2 oo[10 14 ft- 1 M s ] 
r-200 [h- 1 Mpc] 
T spcc [keV] 

N H [10 20 cm- 2 ] 

llcfcrcnccs.-(i;Oirardi ct al. ( 1998);(2;Heipricli & Bohringer (2002);(3;ikcbc et al. (2002)" 
(4 s Shibate"et al.l H2001):(5 s Dickcv & Lockman (1990) 



4 - y '~0.57 

1.64 
8.071°;?,? 

0.93 



hi. 22 
-1.25 

0.89 

o fi n + 0.05 

8.1 



3.82: 



0.17 



6.42" 



0.06 



4.79 



13.9 



2.04 

1.26 

t 9 -+0.04 
z -°-0.05 

2.58 



I.97Z 



1*2 
1*2 
3,U 
5 
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Fig. 1. — Scatter plots of temperatures and densities of SPH particles. Red and blue points indicate particles with unphysical temperat 
and densities which are removed in computing the X-ray emission. 
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Fig. 2. — The examples of the mock spectrum of two simulated clusters. The top panel shows the results for "Virgo" and the bottom 
panel for "Perseus." Black marks are the mock spectra. Green line provides the best-fit spectrum, Red line is that of a single temperature 
thermal model with the temperature T = Tow 11 '™. Each panel has two residuals in terms of sigmas with error bar of size one. Upper one 
is the residual of the best-fit spectrum (T = T ape c). Lower one is that of the thermal model with the emission weighted 

temperature (T = Te™' m ). 
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2 4 6 

T sim ' m [keV] 

ew l j 

Fig. 3. — The upper panel shows the relation of T ew ' to T spC c (asterisk shows IDEAL, square ACIS, and plus MOS) and 7^ lm ' m 
(open circle). The lower panel shows the difference between T apcc and T*S xa ' xa i 8 = 100(T spcc /T a S j m ' m — 1)[%], where T ape c is the best-fit 
temperature of the mock spectra and J' slm ' m ; s given by equation JHJ. 




'/'zoo '/'zoo 

Fig. 4. — The radial profile of simulated clusters. Square provides the (electron number) density profile and dashed line is its fitting line 
assuming the beta model. Plus shows the temperature profile and solid line is its fitting line assuming the polytropic model. Each square 
and plus point corresponds to the shell with a width of 67/i -1 kpc. 
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Fig. 5. — The distribution of Sj- = T/(T) and 8 n = n/{n). Thick solid lines present the distribution throughout the mesh points of within 
r = r200- Dashed lines are fitting lines of the log-normal distribution. Thin solid lines are the distribution of the shells each 67/i _1 kpc 
distance from the center. Each color indicates different radial interval: r < 335?i~ 1 kpc (green), 335/i _1 kpc < r < 670h~ 1 kpc (blue), and 
r > 675h — 1 kpc (cyan). 
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Fig. 6. — The shape of k li (o"ln t) adopting a = 0.75. In the case of ctln y = 0.1 and 0.3, k li ~ 0.99 and 0.89, respectively (See Table 
0. 
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Fig. 7. — The bias due to the radial profile ft RP (eq. [28]) assuming the beta model and two temperature models. We assume that the 
density profile is given by the beta-model. We consider three temperature models: the polytropic model (dotted line provides 7 = 1.1, 
dashed line 7 = 1.2) and the cooling cluster model (solid line). We assume the two case of r c /r2oo- One is r c /r2oo = 1/10 (black line). 
Another is r c /r 200 = 1/40 (red line). 
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Fig. 8. — The emission weighted and spectroscopic like temperature provided our model and the simulation. Dashed line shows 
T Bhn,m = jrmodel or ysim.m = ymodel Do tted lines show T e s 4 m,m /T e ™ odcl - 1 = ±0.1 or T*™'™ /T™° dcl - 1 = ±0.1. In all clusters except 
"Perseus", the temperatures of the model reproduce that of the simulation within 10 percent. 
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Fig. 9. — The bias factor n provided our model and the simulation. Squares show K mo del- Asterisks and crosses show re RP and k li which 
are calculated from our model. In all cases, re mo dcl is kept within ~ 10 %. 



